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Abstract 

The Teukolsky master equation and its associated spin-weighted spheroidal harmonic 

decomposition simplify considerably the study of linear gravitational perturbations of the 
Kerr(-AdS) black hole. However, the formulation of the problem is not complete before we 
assign the physically relevant boundary conditions. We find a set of two Robin boundary con- 
ditions (BCs) that must be imposed on the Teukolsky master variables to get perturbations 
that are asymptotically global AdS, i.e. that asymptotes to the Einstein Static Universe. In 
the context of the AdS/CFT correspondence, these BCs allow a non-zero expectation value 
for the CFT stress-energy tensor while keeping fixed the boundary metric. When the rota- 
tion vanishes, we also find the gauge invariant differential map between the Teukolsky and 
the Kodama-Ishisbashi (Regge- Wheeler— Zerilli) formalisms. One of our Robin BCs maps to 
the scalar sector and the other to the vector sector of the Kodama-Ishisbashi decomposition. 
The Robin BCs on the Teukolsky variables will allow for a quantitative study of instability 
timescales and quasinormal mode spectrum of the Kerr-AdS black hole. As a warm-up for 
this programme, we use the Teukolsky formalism to recover the quasinormal mode spectrum 
of global AdS-Schwarzschild, complementing previous analysis in the literature. 
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1 Introduction and summary 

It is unquestionable that few systems are isolated in Nature and we can learn a lot from studying 
their interactions. Black holes are no exception and the study of their perturbations and inter- 
actions reveals their properties (see e.g. the recent roadmap [I] and review [2] on the subject). 
The simplest deformation we can introduce in a background is a linear perturbation, which of- 
ten encodes interesting physics such as linear stability of the system and its quasinormal mode 
spectrum. Moreover, it also anticipates some non-linear level properties. For example, in the 
collision of two black holes, such as in the coalescence of a binary system, after the inspiral and 
merger phase, the system undergoes a ring down phase where gravitational wave emission is dic- 
tated by the quasinormal mode frequencies. The linear perturbation fingerprints are therefore 
valuable from a theoretical and gravitational- wave detection perspective |li |2j . Perhaps more 
surprisingly, linear analysis of black holes in AdS can be used to infer properties about their 
nonlinear stability [3l[ll[5]. Linear analysis can also infer some properties of (nonlinear) black 
hole collisions and associated gravitational wave emission in the close- limit approximation [6j. 

To study linear gravitational perturbations of a black hole we need to solve the linearized 
Einstein equation. A priori this is a remarkable task involving a coupled system of PDEs. 
Fortunately, for the Kerr(-AdS) black holes (which are Petrov type D backgrounds), Teukolsky 
employed the Newman-Penrose formalism to prove that all the gravitational perturbation in- 
formation is encoded in two decoupled complex Weyl scalars jT] [8] . These are gauge invariant 
quantities with the same number of degrees of freedom as the metric perturbation. Moreover, 
there is a single pair of decoupled master equations governing the perturbations of these Weyl 
scalars (one element of the pair describes spin s = 2 and the other s = —2 modes). In a mode 
by mode analysis, each master equation further separates into a radial and angular equation 
which makes the analysis technically tractable [U [9l \TU[ lllj . (In the absence of rotation, and 



only in this case, we can instead use a similar pair of decoupled master equations for a distinct 
pair of gauge invariant variables proposed by Regge and Wheeler [12j and Zerilli [13] , and later 
recovered and extended by Kodama and Ishibashi |l4j). 

Solving these master equations is not our only task. Like in any PDE system, it is also 
important to assign physically relevant boundary conditions. Without the later, the formulation 
of the problem is not complete. 

In this paper we are interested in linear gravitational perturbations of the Kerr-AdS black 
hole, with a focus on its boundary conditions (BCs). The (extra) motivation to put into a firmer 
basis the linear perturbation problem of the Kerr-AdS system is two- folded. First, the Kerr-AdS 
black hole is known to have linear gravitational instabilities sourced by superradiance [15j-[18j 
and by extremality [19]. Second, in the AdS/CFT duality context, perturbing a (Kerr-AdS) 
black hole in the bulk is dual to perturbing the associated CFT thermal state (with a chemical 
potential) living in bulk boundary. The time evolution of the perturbed black hole maps into 
the time evolution of the thermal state fluctuations and the quasinormal mode spectrum of the 
black hole is dual to the thermalization timescale in the CFT (see e.g |20j-|23]. [2]). 

From the mathematical perspective, the boundary condition choice is arbitrary. We need 
physical input to fix it. Not always but quite often, this leads to a unique choice. We establish 
what are the BCs we need to impose in the Teukolsky master solutions to get perturbations 
that are asymptotically global AdS. To make this statement precise, recall that once we have 
the solution of the Teukolsky pair of master variables (s = ±2), we can reconstruct the metric 
perturbations using the Hertz map [21] -[29]. We get a pair of metric perturbations, one in the 
ingoing radiation gauge (IRG; s = —2) and the other in the outgoing radiation gauge (ORG; 
s = 2). By asymptotically global AdS perturbations we mean that we want the BCs in the 
Teukolsky scalars that yield metric perturbations that decay at asymptotic infinity according 
to the power law found by Henneaux and Teiltelboim |3HI32j. Our task is thus very well defined. 
We have to work out the inverse Hertz map and find how the Henneaux- Teiltelboim metric BCs 
translate into the Teukolsky scalars. 

Before arguing further that this choice should be the physically relevant option, it is illu- 
minating to recall what is the situation in an asymptotically flat system. In this case, the BC 
choice in the Teukolsky scalars amounts to choosing the purely outgoing traveling mode. Intu- 
itively, this is because we are not interested in scattering experiments (where an ingoing mode 
component would be present). Formally, this is because this is the choice that yields a metric 
perturbation preserving the asymptotic flatness of the original Kerr black hole, i.e. conserving 
asymptotically the Poincare group of the Minkowski spacetime. 

A similarly reasoning justifies why the Henneaux- Teiltelboim BCs should be the physically 
relevant boundary condition to be imposed in the Kerr-AdS system [32j. These are the BCs 
that preserve asymptotically the global AdS symmetry group 0(3, 2) and yield finite surface 
integral charges associated with the 0(3,2) generators. Yet an additional reason to single out 
this BC is justified by the AdS/CFT duality. The Kerr-AdS is asymptotically global AdS 
and the CFT lives on the boundary of this bulk spacetime. As desired in this context, the 
Henneaux- Teiltelboim BCs are such that they allow for a non-zero expectation value for the 
CFT stress-energy tensor while keeping fixed the boundary metric. This criterion to select the 
BCs of gauge invariant variables was first emphasized in the context of the Kodama-Ishibashi 
(KI) formalism [14J by Michalogiorgakis and Pufu [33] . They pointed out that previous analysis 
of quasinormal modes of the 4-dimensional global AdS-Schwarzschild black hole using the KI 
master equations were preserving the boundary metric only in the KI vector sector, but not in 
the KI scalar sector of perturbations (here, vector/sector refer to the standard KI classification 
of perturbations). Indeed, previous studies in the literature had been imposing Dirichlet BCs 
on the KI gauge invariant variables. It turns out that in the KI scalar sector keeping the 
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boundary metric fixed requires a Robin BC (which relates the field with its derivate) [33.] • 
Still in the context of AdS/CFT on a sphere, other boundary conditions that might be called 
asymptotically globally AdS (and that promote the boundary graviton to a dynamical field) 
were proposed in [34J. However, they turn out to lead to ghosts (modes with negative kinetic 
energy) and thus make the energy unbounded below |35j . So, the Henneaux-Teiltelboim BCs 
are also the physically relevant BCs for the AdS/CFT where the CFT lives in the Einstein 
Static Universe. 

So, a global AdS geometry with Henneaux-Teitelboim BCs does not deform the boundary 
metric. This is the mathematical statement materializing the pictoric idea that a global AdS 
background behaves like a confining box with a reflecting wall. An interesting observation that 
emerges from our study is that these BCs require that we consider a particular linear combina- 
tion of the Teukolsky IRG and ORG metric contributions. We can interpret this property as 
being a manifestation of the common lore that only a standing wave with a node at the AdS 
boundary can fit inside the confining box. This pictorial notion of a standing wave and node 
is very appealing but, what is the formal definition of a node in the present context? Does it 
mean that we have to impose a Dirichlet BC on the Teukolsky scalars? No. Instead we will find 
the Robin BC (3.9 )-(3.12 ), much like what happens in the scalar sector of the aforementioned 
4-dimensional KI system. An inspection of this Robin BC (pair) immediately convinces us that 
we hardly could guess it without the actual computation. 

At first sight the fact the asymptotically global AdS BC requires a sum of the Teukolsky 
IRG and ORG metric components is rather surprising and, even worrying. Surprising because 
in the asymptotically flat case we just need to use the outgoing contribution. Eventually wor- 
rying because it is known that in Petrov type D backgrounds the two Teukolsky families of 
perturbations (s = ±2) encode the same information, once we use the Starobinsky- Teukolsky 
identities [361 EZl EH 13 ED] that fix the relative normalization between the s = ±2 Weyl scalars 
perturbations. Our result is however not in contradiction with this property. Indeed, the pre- 
vious statement implies that the most general solution of the s = ±2 master equations contain 
the same information but says nothing about the BCs. This just highlights that the differential 
equations and the BCs are two distinct aspects of the problem, which is not a surprise. Once we 
find our BCs, in practical applications, we just need to study the s = 2 (say) Teukolsky sector 
of perturbations. We believe that an infinitesimal rotation of the tetrad basis should allow to 
derive our results using only the outgoing gauge (say) , although at the cost of loosing contact 
with the standing wave picture. 

We have already mentioned that perturbations for static backgrounds (with global AdS 
and global AdS-Schwarzschild being the relevant geometries here) can be studied using the 
Kodama-Ishibashi (KI) gauge invariant formalism [H] (i.e. the Regge- Wheeler— Zerilli formal- 
ism jT2| I13j). On the other hand, the Teukolsky formalism also describes these cases when 
rotation is absent. Therefore, the two formalisms must be related in spite of their differences, 
although this one-to-one map has not been worked out to date. We fill this gap in the literature. 
The difference that stands out the most is that the KI formalism decomposes the gravitational 
perturbations in scalar and vector spherical harmonics while the Teukolsky formalism uses in- 
stead a harmonic decomposition with respect to the spin-weighted spherical harmonics. These 
harmonics are distinct and ultimately responsible for the different routes taken by the two for- 
malisms. However, both the KI spherical harmonics and the spin-weighted spherical harmonics 
can be expressed in terms of the standard scalar spherical harmonic (associated Legendre poly- 
nomials) and their derivatives. These two maps establish the necessary bridge between the 
angular decomposition of the two formalisms. We then need to work out the radial map, which 
follows from the fact that the metric perturbations of the two formalisms must be the same 
modulo gauge transformations. This gauge invariant differential map expresses the KI master 
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variables (for the KI scalar and vector sector) in terms of the s = 2 (say) Teukolsky master 
field and its first radial derivative and is given in (4.15)-(4.16). To have the complete map 
between the KI and Teukolsky (a = 0) formalisms we also need to discuss the relation between 
the asymptotically global AdS KI BCs and the global AdS Teukolsky BCs. This is done in 
(4.18 )-(4.22 ). The fact that our BCs for the Teukolsky variables match the Michalogiorgakis- 
Pufu BCs for the KI variables is a non-trivial check of our computation in the limit a — )• 0. 
Yet this exercise reveals to be more profitable. Indeed, an interesting outcome is that there is 
a Teukolsky solution/BC that maps to a KI scalar mode/BC and a second one that maps to a 
KI vector mode/BC. This is the simplest possible map between the two formalisms that could 
have been predicted, yet still a surprise. 

With our asymptotically global AdS boundary conditions, the Kerr- AdS linear perturbation 
problem is completely formulated and ready to be applied to problems of physical interest. These 
include finding the quasinormal mode spectrum of Kerr- AdS and the dual CFT thermalization 
timescales and studying quantitatively the superradiant instability timescale of the solution. 
This programme is already undergoing and will be presented elsewhere. Neverthless, as a first 
application, we can recover the quasinormal mode spectrum of the global AdS-Schwarzschild 
this time using the Teukolsky approach. As it could not be otherwise, we recover the previous 
results in the literature both for the KI vector sector (firstly studied in \39\ HOl Wi\ ) and for the 
KI scalar sector first obtained in |33jQ Our results and presentation contribute to complement 
these analysis by plotting the spectrum as a function of the horizon radius, and not just a few 
points of the spectrum. Our analysis focus on the parameter space region of r-(_ /L (the horizon 
radius in AdS units) where the spectrum meets the normal modes of AdS and where it varies 
the most. We will not discuss the asymptotic (for large overtone ^41j-^44j and for large harmonic 
[5]) behaviour of the QN mode spectrum. 

The plan of this paper is the following. Section [2] discusses the Kerr- AdS black hole in the 
Chambers-Moss coordinate frame [56] (instead of the original Carter frame |l5]) that simplifies 
considerably our future discussion of the results. We discuss the Teukolsky formalism, the 
associated Starobinsky- Teukolsky identities and the Hertz map in a self-contained exposition 
because they will be fundamental to derive our results. In Section [3] we find the BCs on the 
Teukolsky variables that yields asymptotically global AdS perturbations. Section |4] constructs 
the gauge invariant differential map between the Teukolsky and Kodama-Ishibashi (Regge- 
Wheeler— Zerilli) gauge invariant formalisms. Finally, in Sectionjsjwe study the QNM spectrum 
/ CFT timescales of the global AdS(-Schwarzschild) background. 



2 Gravitational perturbations of the Kerr- AdS black hole 
2.1 Kerr-AdS black hole 

The Kerr-AdS geometry was originally written by Carter in the Boyer-Lindquist coordinate 
system {T, r, |45j . Here, following Chambers and Moss [56], we introduce the new time 
and polar coordinates {t, x} related to the Boyer-Lindquist coordinates {T, 9} by 

t = ET, x = acose, (2.1) 

^As noticed in [33] the analysis done in [39j-|44] in the KI scalar case does not impose asymptotically global AdS 
BCs and thus we will not discuss further the scalar results of these studies. 
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where H is to be defined in ( 2.3 ). In this coordinate system the Kerr-AdS black hole line element 
reads 



Ar 



(j.2 _|_ T^2~) ^2 
„2 I ,,2^ 



dt 



a 

^2 I ,,2\ 



A. 



(r^ _l_ T^2^ ^2 



dt 



0? + 



(2.2) 



where 



A,. = (a^ + r^) (^1 + ^^ -2Mr, A^ = (a 



1--^. (2.3) 



The Chambers- Moss coordinate system {t, r, x, 4>} has the nice property that the line element 
treats the radial r and polar x coordinates at an almost similar footing. One anticipates that this 
property will naturally extend to the radial and angular equations that describe gravitational 
perturbations in the Kerr-AdS background. In this frame, the horizon angular velocity and 
temperature are given by 



2ti 



1 + 



^2 
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(2.4) 



The Kerr-AdS black hole obeys R^,^ = — 3L '^gfj,u, and asymptotically approaches global 
AdS space with radius of curvature L. This asymptotic structure is not manifest in (2.2), one of 



the reasons being that the coordinate frame {t, r, x, </>} rotates at infinity with angular velocity 



-a/{L S). However, if we introduce the coordinate change 
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(2.5) 



we find that as r — )• oo (i.e. R — )• oo), the Kerr-AdS geometry (2.2) reduces to 



ds 



AdS 



l + ^\dT^ + 



-^^ + R' {de' + sin^ed^'), 

1 + X2 



(2.6) 



that we recognize as the line element of global AdS. In other words, the conformal boundary 



of the bulk spacetime is the static Einstein universe Rt x S'^: dsl^^y 



lim 



ds 



AdS 



—dT"^ + dQ^ + sin^ Qd^^. This is the boundary metric where the CFT lives in the context of 
the AdS4/CFT3 correspondence. 

The ADM mass and angular momentum of the black hole are related to the mass M and 
rotation a parameters through Madm = M/S^ and Jadm = Ma/r? , respectively |l8lll9]. The 
horizon angular velocity and temperature that are relevant for the thermodynamic analysis are 
the ones measured with respect to the non-rotating frame at infinity \A'6\ and given in terms 
of ( 2.4 ) by T/i = H Th and r^/j = H + . The event horizon is located at r = r_|_ (the largest 
real root of A,.), and it is a Killing horizon generated by the Killing vector K = Ot + ^hd^- 
Further properties of the Kerr-AdS spacetime are discussed in Appendix A of [37] . 
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2.2 Teukolsky master equations 

The Kerr-AdS geometry is a Petrov type D background and therefore perturbations of this 
geometry can be studied using the Teukolsky formahsm, which uses the Newman-Penrose (NP) 
framework [TllSlTl]. 

The building blocks of this formalism are: 

• the NP null tetrad = {i, n, m, m} (the bar demotes complex conjugation) obeying the 
normalization conditions i ■ n = —1, m-m = l; 

• the NP spin connection -jcab = e^^e^ V^,eau (with ^y^ab =-lacb ); 

• the associated NP spin coeficients defined in terms of ^^ab ^ 

K = -7311, A =7424, i^=7422, 0" = " 7313, " = 5(7124-7344), /3 = 5 (7433 " 7213) , 

/^=7423, /? = -7314, Vr=7421, T = - 7312, 7 = 5(7122 - 7342), C = 5 (7431 " 7211 ) i 

(2.7) 



• the five complex Weyl scalars {Cabcd ^-re the Weyl tensor components in the NP null basis) 

^0 = (^1313 , *1 = C12IZ , ^2 = (^1342 , ^3 = Ci2A2 , ^"4 = C2424 ; (2.8) 

• and the NP directional derivative operators D = £^dfj,, A = n'^d^, 6 = m'^d^, 6 = rn'^d^ . 
The complex conjugate of any complex NP quantity can be obtain through the replacement 
3 o 4. 

The Kerr-AdS background is a Petrov type D spacetime since all Weyl scalars, except ^'2, 
vanish: ^'q = ^'i = ^'3 = \I'4 = and ^'2 = —M{r — ix)~^- Due to the Goldberg-Sachs theorem 
this further implies that K = X = i' = a = 0. In addition, we might want to set e = by 
choosing £ to be tangent to an affinely parametrized null geodesic i^Vfj^ii, = 0. This was the 
original choice of Teukolsky and Press (when studing perturbations of the Kerr black hole) who 
used the the outgoing (ingoing) Kinnersly tetrad that is regular in the past (future) horizon [38j . 
In the Kerr-AdS case we can work with the natural extension of Kinnersly's tetrad to AdS, and 
this was the choice made in [50j. However, here we choose to work with the Chambers-Moss 
null tetrad defined as [46j, 

which is not affinely parametrized (e 7^ 0). The motivation for this choice is two-folded. First, 
the technical analysis of the angular part of the perturbation equations and solutions will be 
much simpler because this Chambers-Moss tetrad explores the almost equal footing treatment 
of the r, X coordinates much more efficiently than Kinnersly's tetrad. Second, to complete 
our analysis later on we will have to discuss how the metric perturbations hab (built out of 
the NP perturbed scalars) transform both under infinitesimal coordinate transformations and 
infinitesimal change of basis. It turns out that if we work in the Chambers-Moss tetrad, the 
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results will be achieved without requiring a change of basis, while the Kinnersly's option would 
demand it. Again, this simplifies our exposition. 

Teukolsky's treatment applies to arbitrary spin s perturbations. Here, we are interested in 
gravitational perturbations so we restrict our discussion to the s = ±2 spins. Let us denote 
the unperturbed NP Weyl scalars by and their perturbations by with i = 1, • • • ,5. The 
important quantities for our discussion are the scalars and 5^4. They are invariant both 
under infinitesimal coordinate transformations and under infinitesimal changes of the NP basis. 
A remarkable property of the Kerr-AdS geometry is that all information on the most general 
[^linear perturbation of the system is encoded in these gauge invariant variables (5^'o and ^^'4. 
That is, the perturbation of the leftover NP variables can be recovered once 5^o and ^^'4 are 
known. The later are the solutions of the Teukolsky master equations. 

For s = 2 perturbations the Teukolsky equation is 



[D - 3e + e - 4p - ^] (A + /i - 47) 

- [5 + 7f-a-3/3-4r] (5 + 7r-4a) -3^'2|5^'o = 47r75, if s = 2, (2.10) 

while s = —2 perturbations are described by the Teukolsky equation 

[A + 37 - 7 + 4^ + /2] (D + 4e - />) 

- [J- f + ^ + 3a + 47r] ((5- r + 4/3) -3^2! -5^4 = 47171, if s = -2. (2.11) 

The explicit form of the source terms 7(s±)5 that vanish in our analysis, can be found in [8j. 
Next we introduce the separation ansatz 



(2.12) 



f 6^4 = {r- ixV e-^"* e}^'^ S^S^hd , for s = -2 , 

\ S^o = {r- e-^-* e™<^ R^^L^r) S^^}^{x) , for s = 2 . 

Also, define the radial {Vn^Vn} and angular {Cn,Cn} differential operators 

£. = a, + ^+n^, Ci = d^-^+n^, (2.13) 
^x ^x ^x 

where the prime ' represents derivative wrt the argument and 

Kr = r.[ma-uj{a^ + r'^)], K^ = '~[ma-u}{a^-x^)]- (2.14) 



With the ansatz (2.12), the Teukosky master equations separate into a pair of equations for the 



radial R^^j^{r) and angular sl^^^{x) functions, 



' vl.ArV, + 6 (fl - iHcr) - a) Rl-^Jir) = , 
£UA^Cr + 6[^+Eux)+x)sl-,^{x) = 0, for s = -2 , 



^Excluding the exceptional perturbations that simply change the mass or angular momentum of the background 
|55| . The Teukolsky formalism does not address these modes. See Appendix^ for a detailed discussion. 
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(V^^^VI + 6 + iH^r) - a) i?il(r) = , 
(£-iA^£t + 6 - x) + a) Sil(x) = , for 
where we introduced the separation constant 



(2.16) 



Some important observations are in order: 

• First note that the radial operators obey Vn = (Vn)* (where * denotes complex conjuga- 
tion) while the angular operators satisfy Cn{x) = —^n{—x)- 

(—2) 

• Consequently, the radial equation for ii^^^ is the complex conjugate of the radial equa- 
tion for -R^/^, but the angular solutions S^^^ are instead related by the symmetry 

(2) (—2) 

S^j/^ix) = Su)tm^~^)- '^^^ later statement implies that the separation constants are 

(—2) (2) 

such that A^^^ = A^^^^ = A with A being real. 

• The eigenfunctions S^^j^{x) are spin-weighted AdS spheroidal harmonics, with positive 
integer (. specifying the number of zeros, £ — max{|m|, (so the smallest £ is £ = |s| = 2). 
The associated eigenvalues A can be computed numerically. They are a function of w, ^, m 
and regularity imposes the constraints that —i<m<l must be an integer and £ > \s\. 

• We have the freedom to choose the normalization of the angular eigenfunctions. A natural 
choice is 

'sllSdx = l. (2.18) 



2.3 Starobinsky-Teukolsky identities 



Suppose we solve the radial and angular equations (2.16), (2.15) for the spin s = ±2. These 



solutions for R^^j^ and S^^^{x), when inserted in (2.12), are not enough to fully determine the 



NP gauge invariant Weyl scalars (5^*0, 5^'4. The reason being that the relative normalization 
between S^o and ^^'4 remains undetermined, and thus our linear perturbation problem is 



yet not solved [38[ ^ [TU]. Given the natural normalization (2.18) chosen for the weighted 



spheroidal harmonics, the completion of the solution for 5^'o, <5^'4 requires that we fix the relative 
normalization between the radial functions R^^i^ and R^£^- This is what the Starobinsky- 
Teukolsky (ST) identities acomplish [361 133 ESI El EHI- A detailed analysis of these identities 
for the Kerr black hole is available in the above original papers or in the seminal textbook of 
Chandrasekhar [TT]. Here, we present these identities for the Kerr- AdS black hole. 

Act with the operator 'D^_^ArVQVQArVl on the Teukolsky equation (2.16) for -R^^^n ^^'^ 



(—2) 

use the equation of motion for -R^^^ - This yields one of the radial ST identities. Similarly, 



to get the second, act with the operator V^iArVoVoArVi on the Teukolsky equation (2.15) 



(—2) (—2) 

for R^ji^, and make use of the equation obeyed by -R^^^ - These radial ST identities for the 

(2) (—2) 

Kerr- AdS background relate -R^/^ to R^jg^, 



vl.ArVlvlArVlR^Jl = C.,rI-^ , 

(-2) _ ^* „(2) y^-^^) 



(jjim ' 

V-iArVoVoArV.Rl^i = Cl,R^,^ , 



where we have chosen the radial ST constants {Cst,C*t} be related by complex conjugation. 
This is possible because, as noted before, the -R^^^ solutions are related by complex conjugation. 
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To get the angular ST identities, act with the operator A^£q/^qA^£| on the Teukolsky 

. . (—2) 

equation (2.16) for S^/^ (and use the equation of motion for S^^J^), and act with C-iA^CoCoAf^Ci 



on the equation (2.15) for sj^^^ (and use the equation for S^^^^). This yields the pair of ST 
identities, 



(2.20) 



(±2) 

Since the equations for S^gJ^ are related by the symmetry x —X: the ST constant /Cst on the 

(±2) 

RHS of these ST identities is real. Moreover, because S^^^^ are both normalized to unity — see 



(2.18) — the ST constant is the same in both angular ST identities. 



2 A h 

To determine |Cst| we act with the operator of the LHS of the first equation of ( 2.19 ) on the 



second equation and evaluate explicitly the resulting 8*^ order differential operator. A similar 



operation on equations (2.20) fixes /Cg^. We find that 

|CstP = ICi + UAM^uj^E"^ , (2.21) 
/C^^ = (A + 2f + SXE'^au [(6 + 5A) (m - aco) + 12auj] + U4E'^a'^uj'^{m - aujf 

+ + 2) (A - 6) + 12H2(m - auj) [2mA - aw (A - 6)]] + • (2-22) 

This fixes completely the real constant /Cst (we choose the positive sign when taking the square 
root of /Cg^ to get, when a — )• 0, the known relation between the s = ±2 spin-weighted spherical 
harmonics) but not the complex constant Cgt- However, we emphasize that to find the asymp- 
totically global AdS boundary conditions in next section, we do not need to know Cst, just 
/Cst- Moreover, we do not need the explicit expression for Cst to construct the map between the 
Kodama-Ishibashi and the a = Teukolsky formalisms of Section |4j 

Neverthless we can say a bit more about the phase of Cst- Recall that in the Kerr case, 
finding the real and imaginary parts of Cst requires a respectful computational effort which 
was undertaken by Chandrasekhar [TU] (also reviewed in sections 82 to 95 of chapter 9 of the 
textbook [TT])- A priori we would need to repeat the computations of [lOJ, this time in the 
AdS background, to find the phase of Cst in the Kerr- AdS background (which was never done 
to date). However, if we had to guess it we would take the natural assumption that Cst is given 



by the solution of (2.21) that reduces to the asymptotically flat partner of [TO] when L — t- oo, 



Cst = Ci + iC2 with Ci = /Cst , C2 = -UMcoE. (2.23) 

However, we emphasize again that this expression must be read with some grain of salt and 
needs a derivation along the lines of to be fully confirmed. 

Having fixed the ST constants we have specified the relative normalization between the 
Teukolsky variables 6^0 and ^^'4. We ask the reader to see Appendix [A] for a further discussion 
of this issue. 



2.4 Metric perturbations: the Hertz potentials 

In the previous subsections we found the solutions of the Teukolsky master equations for the 
gauge invariant Weyl scalars of the Newman-Penrose formalism. We will however need to know 
the perturbations of the metric components, h^i, = 5g^u- These are provided by the Hertz map, 
h^u = h^i,{ipH), which reconstructs the perturbations of the metric tensor from the associated 
scalar Hertz potentials ipn (in a given gauge) [24j-[29j. The later are themselves closely related 
to the NP Weyl scalar perturbations and (5\I'4. 
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In the Kerr-AdS background, the Hertz potentials are defined by the master equations they 
obey to, namely, 



(-2) 



0. 



[(A + 37 - 7 + + 4e + 3p) - (5 + /3 + 3a - r)((5 + 4/3 + 3r) - 3^2] 
[{D-Se + e- p)(A - 47 - 3^) - ((5 - 3/3 - a + 7f)(^ - 4q - 37r) - 3*2] V'^^ = . 

Introducing the ansatz for the Hertz potential 

-2, 



s = 
s = 2. 



(2.24a) 
(2.24b) 



(2.25) 



into (2.24) (in the Kerr-AdS background), we find that and S^^^ are exactly the solutions 



of the radial and angular equations (2.15) and (2.16). This fixes the precise map between the 



Hertz potentials and the NP Weyl scalar perturbations. 

(—2) (2) 

The Hertz map is such that the Hertz potentials -0/^ and V'h generate the metric pertur- 
bations in two different gauges, namely the ingoing (IRG) and the outgoing (ORG) radiation 
gauge, defined by 



IRG : £^/i„ 



0, g'^'-h, 



0. 



ORG : V = 0' 9'"'^ 



0. 



(2.26) 



The Hertz map is finally given b}|^ 

hf^ = {kp^m^) [(^ + 3e + e - p + p) (,5 + 4/3 + 3t) + (J + 3/3 - Q - r - 7f) (D + 4e + 3p)] 
-£^4 ((5 + 3/3 + a - r) ((5 + 4/3 + 3r) - m^m^ (D + 3e - e - p) (Z) + 4e + 3p)}V'{7^^ 
+C.C. , (2.27) 

h^f^ = {^(^m^) [(^ + ^ - 3a + r + vr) (A - 47 - 3p) + (A - 37 - 7 + // - /i) (5 - 4a - 37r)] 

(5 - ^ - 3a + vr) (5 - 4a - 37r) - m^m^ (A - 37 + 7 + p) (A - 47 - 3/i)}V'S^^ 

(2.28) 



-l-c.c. 



We have explicitly checked that (2.27) and (2.28) satisfy the linearized Einstein equation (see 
also footnote [S]) . 

It is important to emphasize that the Hertz map provides the most general metric pertur- 
bation with ^ > 2 of the Kerr-AdS black hole [81 [9l [TOl |27] . We defer a detailed discussion of 
this observation to Appendix [A} 



3 Boundary conditions for global AdS perturbations of Kerr- 
AdS 

We start this section with a brief recap of the Teukolsky system which emphasizes some of its 
properties that are essential to discuss the asymptotic boundary conditions. 

The gravitational Teukolsky equations are described by a set of two families of equations, 
one for spin s = 2 and the other for s = —2. In Petrov type D backgrounds, these two 
families encode the same information, once we use the Starobinsky- Teukolsky identities that fix 

(s) 

the relative normalization between both spin-weighted spheroidal harmonics S^^ and between 



Note that (2.281, whose exphcit derivation can be found in an Appendix of |29| . corrects some typos in the map 



first presented in [25) , 
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both radial functions n^i^- Indeed, modulo the ST relative normalization, the two radial 
functions are simply the complex conjugate of each other, and the two angular functions are 
related by ^^^/^(x) = '^ujem(~^)- This is a consequence of the fact that the Teukolsky operator 
acting on is the adjoint of the one acting on ^^'4. 

The upshot of these observations, with relevance for practical applications, is that the 
Teukosky system in Petrov type D geometries is such that we just need to analyze the s = 2 

sector (for example) to find all the information, except BCs, on the gravitational perturbations 

(2) 

(excluding modes that just shift the mass and angular momentum). In other words, given 

(2) 

^uiim constants we can reconstruct all the s = —2 Teukolsky quantities. 

Were we discussing perturbations of the asymptotically flat Kerr black hole and this section 
on the boundary conditions would end with the following single last observation. Being a sec- 
ond order differential system, the gravitational field has two independent asymptotic solutions, 
namely, the ingoing and outgoing traveling modes. Since we are not interested in scattering 
experiments, the EC (that preserves asymptotic flatness) would be fixed by selecting the purely 
outgoing BC. For practical purposes, we would definitely just need to study the s = 2 Teukolsky 
system of equations. 

The situation is far less trivial when we look into perturbations of Kerr-AdS. This time the 
second order differential system has two independent asymptotic solutions that are power laws 
of the radial variable. The BC to be chosen selects the relative normalization between these 
two solutions. What is the criterion to make this choice? This will be made precise in the 
next subsection. Before such a formal analysis we can however describe it at the heuristic level. 
Basically we want the perturbed background to preserve the asymptotic global AdS character 
of the Kerr-AdS background. Global AdS asymptotic structure means that the system behaves 
as a confining box were the only allowed perturbations are those described by standing waves. 
Standing waves on the other hand can be decomposed as a fine-tuned sum of IRG and ORG 
modes such that we have a node at the asymptotic AdS wall. With this brief argument we 
conclude that to find the asymptotic global AdS BC we necessarily need to use the information 
on both the IRG and ORG Teukolsky metric perturbations, i.e. the BC discussion will require 
using information on both spins. Once we find it, it is still true that the spin s = 2 sector of 
the Teukolsky system encodes the same information as the s = —2 one, and we will be able 
to study the properties of perturbations in Kerr-AdS using only the s = 2 sector (say). (Note 
that an infinitesimal rotation of the tetrad basis should allow to derive our results using only 
the ORG, say). 

So we take the most general gravitational perturbation of the Kerr-AdS black hole to be 
given by the sum of the ingoing and outgoing radiation gauge contributions as written in (2.27) 
and (2.28). (By diffeomorphism invariance, this solution can be written in any other gauge 
through a gauge transformation). The physically relevant perturbations are those that are 
regular at the horizon and asymptotically global AdS. In this section we find one of our most 
fundamental results, namely the BCs we need to impose on our perturbations. 



3.1 Definition of asymptotically global AdS perturbations 

When considering linear perturbations of a background we have in mind two key properties: 
the perturbations should keep the spacetime regular and they should be as generic as pos- 
sible, but without being so violent that they would destroy the asymptotic structure of the 
background. To make this statement quantitative, in the familiar case of an asymptotically 
Minkowski background, the appropriate boundary condition follows from the requirement that 
the perturbations preserve asymptotically the Poincare group of the Minkowski spacetime j30j . 
For the AdS case, Boucher, Gibbons, and Horowitz [3T] and Henneaux and Teiltelboim [32] 
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have defined precisely what are the asymptotic BC we should impose to get perturbations that 
approach at large spacelike distances the global AdS spacetime. The main guideline is that per- 
turbations in a global AdS background must preserve asymptotically the global AdS symmetry 
group 0(3,2), much like perturbations in a flat background must preserve asymptotically the 
Poincare group of the Minkowski spacetime. More concretely, asymptotically global AdS space- 
times are defined by BCs on the gravitational field which obey the following three requirements 

(1) they should contain the asymptotic decay of the Kerr-AdS metric; 

(2) they should be invariant under the global AdS symmetry group 0(3, 2); 

(3) they should make finite the surface integral charges associated with the 0(3, 2) generators. 
If we work in the coordinate system {T, i?, 0,$}, where the line element of global AdS is 



given by (2.6), the metric perturbations that obey the above BCs behave asymptotically as f32]: 



hT^l -- 

hRR 

hR^i -- 



-Ft^{T,Q,^) + 0{R-^) 
1 

1 

' W 
1 

1 

R 



for /i = T, e, ^> , 



FTi?(r,e,$) + o(fi-5) 

FiiR{T,e,1>) + 0{R-^) 

Fr^{t,q,<^>) + o{r-'^) 



F^,(r,e,$) + o(i?-2) 



for /i = 6, $ , 
for /i, = 0, <I> , 



(3.1a) 
(3.1b) 
(3.1c) 
(3.1d) 
(3.1e) 



where Ffj,u{T, 0, <1>) are functions of {T, 0, <!>} only. 

These BCs are defined with respect to a particular coordinate system. Consider a generic 



infinitesimal coordinate transformation 



where ^ is an arbitrary gauge vector field. 



Under this gauge transformation the metric perturbation transforms according to 



(3.2) 



which we can use to translate the BCs (3.1 ) in the {T, R, 0, <!>} frame into any other coordinate 
system, so long as ^ decays sufficiently fast at infinity. 



3.2 Boundary conditions for asymptotically global AdS perturbations 

Modulo gauge transformations, the most general perturbation of linearized Einstein equations 
in the Kerr-AdS background can be written as 



h 



111/ 



ORG 



(3.3) 



where h^^'^ and h^^'^ are determined by the Hertz maps (2.27) and (2.28), with the Hertz po- 



f±2) 

tentials V')^ defined in (2.25) and the associated Teukolsky functions obeying the equations of 
motion (2.16) and (2.15). Note that the relative normalization between these two contributions 



is fixed by the Starobinsky- Teukolsky treatment. 



Solving the radial Teukolsky equations (2.16) and (2.15) at infinity, using a standard Frobe- 

(±2) 

nius analysis, we find that the two independent asymptotic decays for are 



A 



(2) 



L 



+ A 



(2) 



B 



(-2) 



L 



B 



(3.4) 
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where the amphtudes A 



(^) 
± 



(—2) 

} are, at this point, independent arbitrary constants. 



Our task is to find the BC we have to impose in order to get a perturbation hnu that is 



asymptoticaUy global AdS. That is, we must find the constraints, A 



these amplitudes have to obey to get the Henneaux-Teiltelboim decay (3.1). We will find that 



that 



the most tempting condition, where we set to zero the leading order term in the expansion, 



A 



is) 



0, is too naive and does not do the job. Note that it follows from (2.5) that for large R, 

1 1/2 



or r, one has R ^ r [(L^ — x^)/-^^^]^^^ and cos© ~ x [-^^'^'/(-^^ ~ X^)]^^^- Therefore, to get 
the asymptotically global AdS decay of h^^ in the {t, r, x, 0} coordinate system we can simply 
replace {r,i?,G,^>} 



{t,r,x,4>} in (3.1). 



In ■, we can express s[ji^{x) as a function of S]^l^{x) and its derivative using the first 



(2) 



C_2] 

eliminates -S^^^ from (3.3). At this stage, we could also immediately replace the ST constant 
/Cst by its expression (2.22). Instead, we choose to keep it unspecified until a later stage in our 
computation. 

The explicit expression of h^^^ + when we introduce into ( |2.27D and ( |2.28D 

contains order terms but no other higher power of r. Our first task is to use all the gauge 



Starobinsky-Teukolsky (ST) identity in (2.20) and the angular equation of motion (2.16). This 



freedom (3.2) to eliminate, if possible, these 0{r ) terms and all lower power law terms that 



are absent in the asymptotically global AdS decay (3.1). The gauge parameter compatible with 
the background isometrics is ^ = e~*'^*e*™''^^^(r, x)dx'^. A simple inspection of V(fj,^u) concludes 
that the most general components of the gauge vector field, that can contribute up to 0{r'^) 
terms, can be written as the power law expansion in r: 



3=0 



3=0 



3=0 



3=0 



(3.5) 

Inserting this expansion and (3.4) into (3.3), we find that there is a judicious choice of the 



functions C^*''(x) such that we can eliminate most of the radial power law terms that are absent 



in the several metric components of (3.1) (the expressions are long and not illuminating). More 
concretely, we are able to gauge away all desired terms but the 0{r'^) contribution in the 
components /lyy, h^A, and hsd,- 



At this point, having used all the available diffeomorphism (3.2), we find ourselves at a key 
stage of the analysis. To eliminate the undesired leftover O(r^) contributions we will have to fix 



the BCs that the amplitudes introduced in (3.4) have to obey to guarantee that the perturbation 
is asymptotically global AdS. There are two conditions that eliminate simultaneously the ©(r^) 
terms in h^^, h^^ and /i^^. One is the coefficient of a term proportional to 'S'^/m(x); and the other 
is proportional to 9;^5'^£m(x)r Clearly, these two contributions have to vanish independently. 
We can use them to express, for example, the amplitude B]^ and the ST constant Kst in 

(—2) (2) 

terms of the other amplitudes B_ ,Aj^\ perturbation parameters uj,m,X and the rotation 



Note that in this computation we use the angular equation of motion (2.161 to get rid of second and higher 
derivatives of S' 



(2) 
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background parameter a (the mass parameter M is absent in these expressions): 



AamLE(2iA'L^ + 5LE ojA'i^A + 2c? (-6 + A + ^l^'E^u?^ A^J} 



JC, 



B 



(-2) 



St 



.21 r2^2^2 



(3.6) 

|a2(2 + A)2l6 + 8a A(6 + bX)mL^ r?Lo - lUa^mL^r?u: (-2 + A + 2L^r?uj'^) 
A [-12 + (-4 + A)A + 2W] + 2 [(6 - 5A)A + ISm^] L 
36 - 12A + A^ - 48Am2 + 12 [A - 6 (l + rr?)] L^E?uj^ + 36L'^H'^a;^ 
+48a^m2 (2A + 2,L^r?uj'^) | , (3.7) 



where we have defined 



A(2 + X)a'l^ -4{l + X)LEuj( iXA'^^^' + 2LE ujA^_^> ) + 4i(2 + 3A)L^H^lj^A 



,(2) 



(2) 



'3^3, ,3 4(2) 



+8L7hV (a^^^ -iLEooA^^A -AamL'E 



ULEujA'i^ 



+2a^L 



2[2iLE ujA^l^ -A'^^Ai'i + X-^, L^EW) + 3XA^^^ 



(3.8) 



At this stage we finahy introduce the explicit expression for the angular Starobinski-Teukolsky 



constant, namely, /Cst is given by the positive square root of (2.22). In addition, we also use 

(±2) 

the property that the radial -R^^^ solutions are the complex conjugate of each other. In these 



conditions we find that conditions (3.6)-(3.7) are obeyed if and only if 

,(2) 



A' 



(3.9) 



with two possible solutions for rj, that we call r]s and r]v, for reasons that will become clear in 
the next section. These define the BCs we look for. To sum up, the two possible BCs on the 



take the form (3.9) with 



Teukolsky amplitudes, defined in (3.4), that yield an asymptotically global AdS perturbation. 



1) V = Vs - 

2) r] = r]v 



Aq-. 
A2 

Aq + ^/A^ 
Ao 



or 



(3.10) 
(3.11) 



where we have introduced 

Ao = 2a^(A - 6) - 8(A + 1)LWe'^ + 8LVS^ + [A(A + 2) - 4 E'^au [5{m - au) + 2auj]] , 
Ai = 4a^(A - 6)2 + L^a2(A + 2f + 48(A + 6)a^E'^L'^uj{m - au) + 8A(5A + 6)(m - auj)L^E^aio 

+4.a^L^ X [-12 + (A - 4)A + 24(m - aujfE^] + UE^L^oj"^ [2X + 3(m - au:fE^] , 

A2 = 4LH [2am + L^u (2 + A - 2LVh2)] . (3.12) 
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Note that the BCs do not depend on the mass parameter M of the background black hole 
(neither does the ST constant /Cst)- 

The metric of the Kerr-AdS black hole asymptotically approaches that of global AdS. The 
boundary conditions (3.9)-(3.10) are the most fundamental result of our study: perturbations 
obeying these BCs are the ones that preserve the asymptotically global AdS behavior of the 
background. These are also natural BCs in the context of the AdS/CFT correspondence: they 
allow a non-zero expectation value for the CFT stress-energy tensor while keeping fixed the 
boundary metric. The reader interested in different BCs that allow, e.g. for a dynamical 
boundary metric, can start from the respective asymptotic metric decay that replaces (3.1) and 
work out the above procedure to get the associate BCs on the Teukoslsky variables. 



3.3 Horizon boundary conditions 

At the horizon, the BCs must be such that only ingoing modes are allowed. 
A Frobenious analysis at the horizon gives the two independent solutions, 



Kim ~ 



A. 



+ Aout (r - r+) 



+ 0{{r 



(3.13) 



where ^in, Aout are arbitrary amplitudes and Q,h,Th are the angular velocity and temperature 



defined in (2.4). The BC is determined by the requirement that the solution is regular in ingoing 



Eddington-Finkelstein coordinates (appropriate to extend the analysis through the horizon) and 
demanding regularity of the Teukolsky variable in this coordinate system. This requires that 
we set Aout = in (3.13): 



R 



is) 



Ain (r - r+) 



+ 0((r-r+)2) 



(3.14) 



4 Map between Teukolsky and Kodama-Ishibashi formalisms 

(a = 0) 

In the previous section we found the boundary conditions we need to impose on the solution 
of the Teukolsky master equation to get gravitational perturbations of the Kerr-AdS black 
hole that preserve the asymptotically global AdS behavior of the background. The Kerr-AdS 
family includes the global AdS-Schwarzchild black hole and the global AdS geometry as special 
elements when we set, respectively, a = and a = = M. Thus, our BCs also apply to 
perturbations of these static backgrounds. 

On the other hand, perturbations of the global AdS(-Schwarzchild) backgrounds were al- 
ready studied in great detail in the literature using the Kodama-Ishibashi (KI) gauge invariant 
formalism. In four dimensions, the KI formalism reduces exactly to the analysis firstly done by 
Zerilli and Regge and Wheeler (in the L — >■ oo case). Indeed, the KI vector master equation is 
the Regge- Wheeler master equation for odd (also called axial) perturbations |12j . and the KI 
scalar master equation is the Zerilli master equation for even (also called polar) perturbations 

m ■ 

Clearly, there must be a one-to-one map between the Kodama-Ishibashi and the Teukolsky 
formalisms (when a = 0). This map was never worked out so we take the opportunity to find 
it[^ Actually, this task will reveal to be quite fruitful since we will find some remarkably simple 
connections. 

^For an earlier discussion of the relation between the Teukolsky and the Regge- Wheeler— Zerilli variables in the 
asymptotically flat case see |51| . 
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Recall the main difference between the Kodama-Ishibashi and the Teukolsky formalisms. 
The former is well suited only for backgrounds that are locally the product of a 2-dimensional 
orbit spacetime (parametrized by the time and radial coordinates) and a base space fC. If this 
is the case, we can do a harmonic decomposition of the perturbations hab according to how they 
transform under coordinate transformations on fC. This is certainly the case of the global AdS(- 
Schwarzchild) backgrounds where the base space is a sphere, K, = S"^, and we can introduce a 
sp/ierica/ harmonic decomposition of gravitational perturbations. (Unfortunately, the Kerr-AdS 
geometry cannot be written as a local product of two such spaces and the KI formalism does 
not apply to it). On the other hand, the Teukolsky formalism uses a harmonic decomposition 

is) 

with respect to the spin-weighted spherical harmonics S*^/^ • 

These two harmonic decompositions are distinct and responsible for the differences between 
the Teukolsky and KI formalisms. We can however write uniquely the scalar /vector KI har- 
monics in terms of standard scalar spherical harmonics (associated Legendre polynomials), and 
there is another unique differential map that generates the spin-weighted spherical harmonics 
also from standard spherical harmonics. This provides the necessary bridge between the two 
formalisms that leads to their unique map. To appreciate this we will discuss these harmonics 
in detail. Also, to make the KI discussion self-contained, we will briefly review the KI formalism 
in the next subsection, before constructing the desired map in subsection |4.2[ 

So far we have not discussed the role of the boundary conditions (BCs) in this map. We 
found a set of two distinct BCs for the Teukolsky solution. Quite interestingly, we will see 



that Teukolsky perturbations with EC (3.10) maps to the KI scalar modes while Teukolsky 



perturbations with BC (3.11) generates the KI vector modes. 



4.1 Kodama-Ishibashi gauge invariant formahsm 

In the Kodama-Ishibashi (KI) formalism [14J, the most general perturbation of the global AdS(- 
Schwarzchild) geometries is decomposed into a superposition of two classes of modes: scalar and 
vector. Scalar and vector modes are expanded in terms of the scalar S(x, (/>) and vector Vj(x, (j)) 
harmonics that we review next (x = cos^; see (|2.1[)). 



Scalar perturbations are given by [13] 

hl^ = fab^, hl, = rfaSi, h|J=2r^HL^^JS + HTS,j), (4.1) 

where (a, 6) are components in the orbit spacetime parametrized by {t,r}, are legs on 

the sphere, {fab, fa, Ht, Hi} are functions of (i,r), and S is the KI scalar harmonic, Sj = 
—As ^^'^DiS, Sij = DiDjS-\-^'yijS, jjk is the unit radius metric on S'^ and Dj is the associated 
covariant derivative. Assuming the ansatz S{x,(j)) = e™'^Y^™(x), the scalar harmonic equation 
(A52 + As) S = (A52 = -f^^DjDk) reduces to 

a. [(1 - x^) d.Yi^x)] + (as - ^) y/:(x) = . (4.2) 

Its regular solutions, with normalization J^^ d<p f^^dx |S|^ = 1, are 




(4 + 1), 4 = 0,1,2,---, |m|<4. 

where P^{x) is the associated Legendre polynomial. Hence, the KI scalar harmonic S(x,(/)) is 
the standard scalar spherical harmonic Y^{x, (p). 
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On the other hand, the KI vector perturbations are given by [U] 

Kb = 0, Ki = rhaW,, /i-- = -2A-V2/i^Z)(,V,) , (4.4) 
where {/la, /it} are functions of {t,r) and the KI vector harmonics Vj are the solutions of 

(A52 + Av) = , DjY^ = (Transverse condition). (4.5) 
The regular vector harmonics can be written in terms of the spherical harmonic Y^{x, (j)) as 

V,- dx^ = -^^y™(x, cP) dx+{l- x^) d,Y,^{x, 0) # , (4.6) 
with Av = ^v(^v + l) - 1, U = l,2,---, |m|<^v 

With this harmonic decomposition, the linearized Einstein equation reduces to a set of two 

(7) 

decoupled gauge invariant KI master equations for the KI master fields ^/ which can be written 
in a compact form as, 

U\ f\ 2M 
^2 - ^] <^y'{t,r) = 0, where f = l + — , and j = {s,v}, (4.7) 



f J ^ ^ ' r 

for scalar (s), and vector (v) perturbations. Here, ^2 is the d'Alembertian operator in the 
2-dimensional orbit spacetime, and the expression for the potentials {Us,Uv} can be found 
in equations (3.2)-(3.8) and (5.15) of KI [T3], respectively. They depend on the properties 
of the background, namely on the mass parameter M and cosmological length L, and on the 
eigenvalues As,Av of the associated (regular) vector, and scalar harmonics defined above. Since 
the background is time-translation invariant, the fields can be further Fourier decomposed in 
time as <^^f\t,r) = e'^'^J^^^^l (r). 

We will need to express KI master fields in terms of the metric functions. Going through 
for u ^ 0, one finds that the gauge invariant KI scalar master field is given by 



[^2^2 + 2/2 - / (A, + rf')] Ht + rf (2/ + rf) + r'fH^ 



(s) ^ -g; r (/,, + A/Hl) + f {2iftr + a; r/,,) 2i [rffl+{rW + f-rff')ft] 
ujf{Xs-2f + rf') ^,u;fi\s-2f + rf') 

2r 

+ As/(As-2/ + r/') 

VA;(As-2/ + r/0' 
On the other hand the gauge invariant KI vector master field is 

= -Jf- (^.hr + rh'r) . (4.9) 



\v ' 

The KI master variables have the asymptotic expansion, 

^Iii,L^o~^o + ^i7 + --- ^^^^^ i = {s,v}. (4.10) 

The linear differential map h^Jf^ = h^Jjj (^^^^^) that reconstructs the metric perturbations (in 
a given gauge) can be read from (4.1) and (|4.8[) (scalar case) and from (4.4) and (4.9) (vector 
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case), if we follow |14| . The requirement that these metric perturbations are asymptotically 
global AdS in the sense described in Section [3] imposes the conditions |33t S] : 



Scalar BC: 
Vector BC: 



^1 -- 



6M 



0. 



(4.11) 
(4.12) 



We can of course consider different asymptotic BCs. For example, past studies on (quasi)normal 
modes have considered the BC <I>o = for scalar modes, instead of (4.11). However, these BCs 



are not asymptotically global AdS, i.e. they do not preserve the boundary metric, as first ob- 
served in [33] . We thus we do not consider them. Other BCs that might be called asymptotically 
globally AdS were studied in [33], but turn out to lead to ghosts (modes with negative kinetic 
energy) and thus make the energy unbounded below |35| . 



4.2 Teukolsky vs Kodama-Ishibashi (Regge- Wheeler— Zerilli) 

Equations (4.3) and (4.6) express the KI scalar S{x,(p) and vector V(x,(/ 



harmonics 

function of the scalar spherical harmonic Y^{x, (p) defined in (4.3). The spin-wheighted spherical 



harmonics S^^{x, (/)) used in the Teukolsky harmonic decomposition are also related to the scalar 
spherical harmonic YP'{x,(j)) through the differential map 




(1 - x^) d. 



x'^ dxYim{x,(j)) 



1 



^S^^{x,c 



Vl-x 
i 



x^ 



x^ 



(4.13) 



and 



These (regular) harmonics obey the angular Teukolsky equations (2.15)-(2.16) with a 
X = £{£ + !) -2. 

The two angular maps just described provide the starting point to bridge the KI and Teukol- 
sky formalisms. We now need the radial map. The metric perturbations of the two formalisms 
must be the same (modulo gauge transformations), i.e. 



-+- 



h 



IRG 

IjLV 



+ h 



ORG 



(4.14) 



where recall that the KI metric on the LHS is given by (4.1), (4.3), (4.4) and (4.6). On the 



other hand, the Teukolsky metric on the RHS is given by (2.27), (2.28), (2.25) and (4.13) 



Fix the LHS of (4.14 ) to have fixed values of KI quantum numbers and £v Then, the most 



natural expectation is that such a KI perturbation is described by a (possibly infinite) sum, in 
the quantum number i, of Teukolsky harmonics (the background is spherically symmetric so 
we can set wlog m = in our discussion; see below). In fact, a mode by mode analysis (using 
properties of internal products) reveals that 1^ = 1^ = 1 (with no sum involved). This simplifies 
considerably the construction of the map. 

Take (4.14) with the identification = = £ with integer ^ > 2. The later inequality 
requires a discussion before proceeding. The KI formalism describes all scalar modes with 
integer -^s ^ (where £s = are perturbations that just shift the mass of the solution and 
£s = 1 is a pure gauge mode), and all vector modes with integer -^v ^ 1 (where 1^ = 1 are 
perturbations that generate just a shift in the angular momentum of the solution) |14l I54j . 
However, the Teukosky quantum number is constrained to be an integer i > \s\ = 2, so the 
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Teukolsky formalism is blind to the modes that generate deformations in the mass and angular 
momentum of the geometry |27l [55] . The simplest way to confirm this is to note that the map 
(4.13) would be trivial for the = 0, 1" modes. 



In these conditions, each metric component on the LHS of (4.14) is proportional to the 



spherical harmonic Y^{x), or to its first derivative, or to a linear combination of both such 
contributions. The same applies to the RHS of (4.14). Matching all the coefficients of these 



angular contributions we can find the radial KI metric functions {fab, fa, Ht, Hl}, {ha, /ir} 



that describe the Teukolsky perturbations. Finally, inserting these radial KI functions into (4.8) 



and (4.9) we express the KI variables ^^£'\r) as a function of the radial Teukolsky functions 



-^i^i'm(^) ^^'^ their first derivative. Note that the above KI metric functions are in a particular 



gauge induced by the Hertz construction. However, (4.8) and (4.9), and thus the associated 



final map, are gauge invariant. To write ^^^'/^r) as a function of a single Teukolsky function, 



e.g. Rl^lm^f), we further use the radial Starobinsky- Teukolsky identity ( |2.19[ ) (which, recall, 
expresses 

We finally end up with the desired gauge invariant map between the KI 



(-2), 



(2'] 

'uiim v'') ^ function of -R^/^ (r) and its first derivative) , and the equation of motion 
(-2), 

'(jjlm ^ 



(2.16) for R 



scalar and vector master fields $ 
1 7(7+2)! 



(s,v) 



(2) 

and the Teukolsky radial function -R^/„ (and its derivative): 
r2 /.^ „ (£-2)! , (£ + 2)1 



/CstCst 7(1^ 2/(Ar + 6M) V"* '*(£ + 2)' 
'■f[{2iu + f') (Ar + 6M)-2A/] (r[ 



+ 



2)! 



+ 12 i Mio 



X < r 



+ 



,(2) 



4f'^ - 2/2 (A + 4 - 2iruj + rf) + (A + 2 + rf) [f^ + Sicof - 2co^) 



(v) 



+f [(A + 2 - 2iujrf + 4ruj{2ru} + i) - r/' (A - 2 + Siroj + 2rf')] 

{£-2)1 {£ + 2)1 



(2) 



(4.15) 



i ^/{£ + 2)\ r 



/CstCst 7(7^ 4/ V ""''{£ + 2)\ {£- 2)\ 
X I [/ (A + 2 - Uru - 2rf') + (/'^ + ?,iojf' - 2^^)] 

+rf {rf -2f + 2icor) (u'^^i^' 



+ 12 i Muj 



(2) 



(4.16) 



In these expressions, A = £{£ + 1) — 2 is the spin-weighted eigenvalue for a = and M is the 
mass parameter of the black hole. Note that for o = 0, the case we are discussing in this 
and following sections, the background spacetime is spherically symmetric. Consequently, the 
radial Teukolsky equations and solutions are independent of the azimuthal quantum number 



(2) 

m. Therefore, henceforth we droped the associated subscript in R p in the maps (4.15)-(4.16) 



(and henceforward). The Starobinsky- Teukolsky angular and radial constants /Cst and Cst are 
given by ( |2.2lD -( |2.23D , which for a = boil down to0 

i + 2)\ 



St 



2)! ' 



--St 



12iMa;. 



(4.17) 



To check our matching we explicitly verify that our KI master fields obey the KI master equa- 
tions (4.7) when Rlj^{r) satisfies the radial Teukolsky equation (2.16), with A = £{£ + 1) 
Note that the map (4.15)-(4.16) is valid also in the asymptotically flat limit L — t- oo. 



2. 



"We ask the reader to revisit the discussion leading to (2.231. For that reason, and without any loss except 



compactness, we prefer to present our forthcoming results without explicitly giving the expression for Cat- 
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To have the complete map between the Teukolsky and KI formaUsm we still need to discuss 



the relation between the asymptotically global AdS KI BCs (4.10)-(4.12) and the global AdS 



Teukolsky BCs (3.9)-(3.11 ). The later, for a = 0, simply reduce to 

,(2) 



1) A 



(2) 



-i 7] A , 



A 



(2) 



-i 7] A, 



(2) 



-Lw 1 + 



A 



A - 2 (L2l^2 _ 1^ 



A 



2Luj 



L(jj . 



It follows from (3.4) and (3.9) that asymptotically 

L 1 



R^'ll ^A 



(2) 



L 



1 



+ 



(4.18) 
(4.19) 

(4.20) 



Consider first the scalar case described by (4.15). Choose the EC to be such that r] = rj. 



as defined in (4.18). In these conditions, inserting (4.20) into the scalar map (4.15) and taking 
its asyn 
$0 + ^1 



its asymptotic expansion we find that it reduces exactly to the KI expression (4.10), ^^^^ 

L 



with 



$1 



6M 



£{e + i) 

,(2) 



(4.21) 



(^+2)! 
(^-2)! 



2Cst y(7T2)! ^(^+1) 



UiMu 



(^ + 2) 
{i-2) 



+ UiMuj 



which matches the global AdS KI EC ( 4.11| ) for scalar modes once we use As = t{£ + 1). We 
see this as one of the most non-trivial tests of our calculations. 

Next, take the vector case described by (4.16). This time select the EC r/ = r/v defined in 



KI expression (4.10) with 



(4.19). Plug (4.20) into the vector map (4.16) and take its asymptotic expansion. We get the 



^•0 = 0, 



SujCst 7(£ + 2)! 



+ mMuj 



Cst + 



{i + 2y. 

{£-2)1 



UiMio 



(4.22) 



which is the global AdS KI EC (4.12) for vector modes. 



So one of the Teukolsky Robin ECs selects the scalar sector of KI perturbations and the 
other selects the vector KI sector. This is the simplest map we could have predicted! 

To sum this section, we found the differential map between the Teukolsky/KI variables and 



ECs. For scalar modes, the differential map is given by (4.15) with the global AdS Teukolsky EC 



(4.18) mapping into the scalar KI EC (4.11) via (4.21). For vector modes, the differential map 



is instead (4.16), and the global AdS Teukolsky EC (4.19) maps into the vector KI EC (4.12) 



through (4.22). Ey continuity, when we turn on the rotation a, we can say that the Teukolsky 



EC (3. 9), (3. 10) generates the "rotating scalar modes", while the EC (3. 9), (3. 11) selects the 



"rotating vector modes". 



5 Global AdS (quasi) normal modes 



The normal modes of global AdS can be studied using the Teukolsky equations and the boundary 



conditions (3.9 )-(3.11 ). We take the opportunity to find these normal mode frequencies since 



the scalar modes are not explicitly derived in the literature. We also revisit, from a Teukolsky 
perspective, the quasinormal mode spectrum of the global AdS-Schwarzschild black hole. 
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Introducing the differential operator definitions (2.13), the s = +2 the Teukolsky equations 
(|5.1l) read 



(2) 



^ ^ (jjim 



+ 



+ 



S. 



(2) 



0, 



+ 



^ - 4ii^; + a;' 



R. 



(2) 



(5.1) 



It follows from (2.1) that in the static case it is appropriate to work in the coordinate system 
{i, r, X, (/)} where x = O'X. In these conditions equations ( |2.16 ) with a = reduces to 

(m + 2x)2" 



dx 



a, (A,,a,i?S^(r)) + 



+ 



:(2) 



iAi' 



Ar 



9r^ 

+ 2(1 + -^)+ Mroj - A 



(5.2) 



The spin weighted spherical harmonic is independent of the mass parameter M and cosmological 
radius L, and can be found analytically. Next, we discuss in detail the regularity analysis that 
leads to the solution (4.13). 

The regular solution at the north pole x = 1 is {2F1 is the standard Hypergeometric function) 



S. 



.('^ -il-x) 



m + 2 
2 



m-2 



with a 



\m 



{1+X) ^ 2^1 

21 + |m + 2| + 1 



1 



1 



d + b) ,\m + 2\ + 1, 



1 



\/4A + 9 . 



(5.3) 



This solution diverges at the south pole x = — 1 as a positive power of (1 + x)~^ (or as ln(l + x) 
in the special case of m = 2) unless we quantize the angular eigenvalue and quantum numbers 
as 

A = ^(£+l)-2, with £ = 2,3,4, \m\ < i , (5.4) 

where we have introduced the quantum number £ with i — max{|m|, |s| =2} giving the number 
of zeros of the eigenfunction along the polar direction. The regular spin s = 2 spherical harmonic 
that solves the angular equation (5.2) is finally 

/0^ l"^ + 2| Im — 2 

41 = (1-^)^(1 + ^)^ 

X 2F1 (^(-2£ + \m-2\ + \m + 2|), ^(2 + 21 + \m - 2\ + \m + 2|), \m + 2\ + 1, 



(5.5) 



with the quantum numbers i,m constrained by the conditions (5.4). This harmonic is valid 
both for the global AdS-Schwarzchild and global AdS backgrounds since the angular equation is 
independent of M. Using the relation between the Hypergeometric function and the Associated 
Legendre polynomial we can rewrite (5.5) as (4.13). 

To study the (quasi)normal modes of these backgrounds we now need to study the radial 
equation (5.2). Since this equation depends on the mass parameter M we need to study the 
cases M = and M > separately. In the next subsection we first find the normal modes 
of global AdS and we study the quasinormal modes of global AdS-Schwarzchild in subsection 
(IS^l). 
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5.1 Normal modes of global AdS 



In the global AdS background (M = 0) 
solution. The solution that is regular at the origin 



the radial Teukolsky equation (5.2) has an exact 
= 0) is 



(2) 



An 1 



I r 

T 



|(Lw-2) 



1 + 



I r 

T 



id 



2FI 



1,^+1 + Lw,2(^ + 1), 



2r 



iL 
(5.6) 



where Aq is an arbitrary amplitude. Asymptotically this solution behaves as 



R. 



(2) 



-Ane 



-Utt{Lui-3£) 



i-2Fi{l 



1 + Lw,2(^+ 1),2) 



+ ^ (Lw + £- 1)2^1 ( 



+ 



{I - l){Luj + I + I 
J+l 



1,£ + 1 + La;,2(^ + 1),2) 
■2Fi(^,^ + 2 + La;,2£ + 3,2) 



+ 



I? 



(5.7) 



(2) 

Comparing this decay with (3.4) we can read the expressions for two amplitudes Aj^ and 

(2) 

A_ . These amplitudes are a priori independent but the requirement that the perturbation is 
asymptotically global AdS constrains them to be related by the BCs (4.19). These BCs quantize 
the frequencies of the perturbations that can fit in the global AdS box, respectively, as 

1) Scalar normal modes of global AdS: ujL = 1 + i + 2p , (5-8) 

2) Vector normal modes of global AdS: ojL = 2 + i + 2p , (5-9) 

where the non-negative integer p is the radial overtone that gives the number of nodes along 
the radial direction and recall that ^ > 2 is an integer. The frequencies (5.8) and (5.9) describe, 
respectively, the scalar and vector normal mode frequencies of global AdS. Note that we are 
adopting the standard KI classification of the scalar/vector perturbations, in the sequence of 
our conclusions of the previous section. Without any surprise, the frequencies (5.8) and (5.9) 
precisely agree with the normal mode frequencies of global AdS that we obtain when we solve 
the KI master equation (4.7) subject to the BCs (4.11) , i.e. <I>i = in the scalar case [33111] 
and $0 = in the vector case [HI E] . 



5.2 Quasinormal modes of global AdS-Schwarzchild 

In this subsection we study some properties of the gravitational quasinormal mode (QNM) 
spectrum of the global AdS-Schwarzschild black hole (GAdSBH). In the AdS4/CFT3 duality, 
this spectrum is dual to the thermalization timescales of the perturbed thermal states of the 
CFT3 living on the sphere, as discussed in [33] (following the detailed analysis of the AdS5/CFT4 
case presented in [23]). We focus our attention in the low- lying QNMs (small radial overtone 
p and harmonic £) because they are expected to dominate the late-time behavior of the time 
evolution. 

Many properties of this gravitational QNM spectrum were already studied with some detail 
in the past. The low- lying KI vector QNMs with global AdS BCs were discussed in |39 1 140 | |33] . 
The asymptotic behavior of these vector modes for large overtone were further analyzed in 
[39]- [S] (see footnote [l]) . On the other hand, the low- lying KI scalar QNMs with global AdS 
BCs were studied in [53] (see also [13]). Finally, the asymptotic behavior of the vector /scalar 
QNMs for large harmonic i was found in a WKB analysis in [5|. 

Our results agree with the vector results of Cardoso-Lemos [5S] and with the results of 
Michalogiorgakis-Pufu [33]. Our conclusions and presentation contribute to complement these 
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previous analysis mainly by plotting the QNM spectrum as a function of the horizon radius, 
and not just a few points of the spectrum. Our discussion will always focus on the parameter 
space region of r_|_/L where the relevant physics is and/or where the spectrum varies the most. 
Given the t — (j) symmetry of the GAdSBH, the QNM frequencies always come in trivial pairs 
of {oj, —a;*}. We just plot the element of the pair with positive real frequency. 



To find the QNM spectrum we solve the Teukolsky radial equation (5.2) numerically subject 



to the asymptotically global AdS BCs, namely, (4.18) in the scalar case and (4.19) in the 



vector case. We use spectral methods to solve the numerical problem, which uses a Chebyshev 
discretization of the grid. We work with the compact radial coordinate < y < 1 and with the 
new radial function qj defined as 



Vector: y = w 1 



(-2) 



r+ (-2) , X 1-^^ f (1 -WttTh r+/ L) , 
Scalar: y = 1 - f , rI/^ = (1 - y)y' ( 1 - ^J^T^^^ " 



(5.10) 



iujL 
2 



The horizon BC (3.14) translates to a simple Neumann BC and the asymptotic BC (4.19) 



yields a Robin BC, in the vector case. For the scalar case, both the horizon BC (3.14) and the 



asymptotic BC (4.18) translate to a Robin BC relating and its derivative. In both cases, we 



get a generalized quadratic eingenvalue problem in the (complex) QNM frequencies uj. We give 
the harmonic i and run the code for several dimensionless horizon radius r+ / L. 

To discuss the results consider first the vector QNM spectrum. We will either plot the 
imaginary part of the dimensionless frequency Im(u;L) as a function of the real part Re(a;L) or 
the QNM real/imaginary parts as a function of the horizon radius in AdS units (r+/L). 




Figure 1: Left panel: Hydrodynamic QNM (which ends in the black point) and the first four 
microscopic QNM curves (that start at the vector normal modes of AdS pinpointed 
as red dots) of the £ = 2 harmonic of the vector QNM spectrum of GAdSBH. Right 
panel: Imaginary part of the hydrodynamic vector QNM as a function of the horizon 
radius in AdS units. For large r^jL^ the data approaches the black curve which 
is the analytical prediction (5.11). The green dots in the hydrodynamic and in the 
lowest-lying (p = 0) microscopic QNM curves are exactly the values taken from Table 
2 of [39]. See text for detailed discussion of these plots. 
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Figure 2: Lowest-lying {p = 0) microscopic vector QNM of the 1 = 2 harmonic of GAdSBH. 

The green dots are exactly the values taken from Table 2 of [39| . The blue (magenta) 
dots are obtained solving numerically the Teukolsky (KI) equations. (See text for 
detailed discussion of these plots) . 



On the Left panel of Fig. [T] we plot the first five low-lying QNMs of the 1 = 2 vector 
harmonic. The points in the vertical line have pure imaginary frequencies and, as r^/L grows 
large, it approaches the black point (0, 0). They are the vector hydrodynamic modes since, in the 
limit r_|_ ^ L, they can be found solving the perturbed Navier-Stokes equation that describes 
the hydrodynamic regime of the CFT3 on the sphere (the associated plasma is conformal, hence 
it has zero bulk viscosity and shear to entropy density ratio rj/s = 1/{Att)) [33l[23]. To leading 
order in the inverse of the horizon radius, this hydrodynamic computation yields the frequency 
I33i 



ujL\ 



+ 21-1 L 



hydro 



+ 



L2 

rl 



(5.11) 



Recall that the hydrodynamic regime requires that the perturbation wavelength is much larger 
than the thermal scale (the inverse of the temperature) of the theory. So this regime is achieved 
when r+/L (and thus the temperature ThL) grows without bound and the perturbation fre- 
quency becomes arbitrary small. This is indeed what happens as we approach the black point 
moving from bottom to top along the vertical line of Fig. [TJ On the Right panel of this fig- 
ure we plot the (imaginary) frequency of the hydrodynamic mode as a function of rj^/L for 



< r^/L < 1000. The black curve describes the analytical hydrodynamic curve (5.11). As pre- 



dicted in [33], we confirm that as rj^/L grows the black curve indeed approaches the numerical 
data. 

Returning to the Left panel of Fig. [T| the red points describe the first 4 radial overtones 
{p = 0, 1, 2, 3) of the i = 2 vector normal modes of AdS; see (5.9). The associated four vector 



QNMs in the Left panel of Fig. [T] are microscopic modes (as oppose to hydrodynamic) because 
as these curves move away from the red points, i.e. as T^jL (and TjjL) grows so does Im(a;L). 
So we never reach the hydrodynamic regime oj <^ Tjj and we need the microscopic theory to 
describe them. This tower of overtones which are continuously connected to the normal modes 
of AdS are often said to form the main series or main sequence of the vector QNM spectrum. 
In the plot we also pinpoint with green dots the points with fixed r+/L = 0.5 and r_|_/L = 1.0 
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in each overtone curve (for the p = and microscopic mode curves we also identify the green 
mode with r+/L = 2). The green dots in the hydrodynamic and in the lowest-lying (p = 0) 
microscopic QNM curves are also exactly the values taken from Table 2 of [;39j • Modes with fixed 
r+/L in the main sequence and for p > 1 scale linearly with r+ / L. To illustrate this property we 
connect with auxiliary dashed black straight line the two set of modes with r+/L = 0.5, 1 that 
we have singled out. The points of the lowest-lying {p = 0) QNM curve of the main sequence 
do not however fit in these lines as is visible in the plot, but apart from this curiosity, the p = 
overtone curve is similar to the higher overtone curves. 





Figure 3: Lowest-lying (p = 0) microscopic vector QNMs of the first 8 harmonics of the 
GAdSBH. From bottom to top we have: £ = 2, 3, • • • ,10. The red dots give the 



vector normal mode frequencies (5.9). The green dots are exactly taken from Table 

2 of m\. 



Further properties of this lowest-lying {p = 0) microscopic QNM curve are displayed in Fig. 
[2J The Left (Right) panel plots the real (imaginary) part of the QNM frequency as a function of 
the horizon radius for r^/L < 2. The curve starts at the vector AdS normal mode frequency (red 
point) and the frequency stays very close to the real axis for rj^/L < 0.2. The real (imaginary) 
part keeps (increasing) decreasing monotonically for larger values of tj^/L (e.g. for r_|_ = lOOL 
one has ojL ~ 266.384 - 184.959 z [39j). The same 3 green points for r+/L = 0.5, 1, 2 of Fig. [l] 
are also plotted here. As an extra check of our numerics, in addition to the blue points, which 



are obtained solving numerically the Teukolsky system with BC (4.19), we also display with 



magenta points the numerical solution of the KI system with BC (4.12). In the sequence of the 



map constructed in Section [4] these two curves have to overlap. This is indeed the case (the 
magenta points are drawn larger to be visible in the plot). 

To analyze the evolution of the vector QNM spectrum as the vector harmonic changes, in 
Fig. [3] we plot the same information as in Fig. [2| and in addition the lowest-lying {p = 0) 
microscopic QNM of the next seven vector harmonics i. More concretely, from bottom to top 
we have the harmonics £ = 2, 3, • • • ,10. From the Left panel, we conclude that there are regions 
in the parameter space (i.e. windows in the range of r_|_/L) where the real part of the frequency 
spectrum is in a first approximation isospectral (i.e. difference between consecutive harmonics at 
constant r_|_ / L is approximately constant) but there are also others where the lowest harmonics 
£ = 2, 3 spoil this property. From the Right panel, we see that the window of r^/L around 
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the global AdS case (r+/L = 0; red dot (0,0)) where the imaginary part of the spectrum is 
approximately flat increases as the harmonic ^ grows. The green dots in the ^ = 2,3 curves are 
exactly the values taken from Table 2 of [39) . 




Re(w L) 

Figure 4: Hydrodynamic QNM (which ends at the black point) and the first three microscopic 
QNM curves (that start at the vector normal modes of AdS pinpointed as red dots) 
of the 1 = 2 harmonic of the scalar QNM spectrum of GAdSBH. The green dots have 
exactly the values taken from Table 1 and 2 of |33j (we have added the r+ = 0.5 green 
points in the p>l curves). (See text for detailed discussion of this plot). 



Consider now the scalar QNM spectrum of the GAdSBH. In Fig. [4] we plot the first four 
low-lying QNMs of the ^ = 2 scalar harmonic. In this case, the spectrum has no pure imaginary 
frequencies, and all the QNM curves are continuously connected to the scalar normal modes 



frequencies (5.8) of AdS here coloured as red dots for the overtones p = 0, 1,2,3. Clearly, the 
p = overtone curve on the left is a special curve. Indeed, in this case as r+/L (and ThL) grows 
the imaginary part of the frequency first decreases but then it has a minimum after which it 
approaches zero (see the black point). That is, we approach the hydrodynamic regime uj <C Th- 
This is thus the scalar hydrodynamic QNM. For the p > 1 the imaginary part of the frequency 
decreases monotonically as r_|_/L: these are microscopic scalar QNMs. Solving the linearized 
hydrodynamic equations on Rt x S"^ for a conformal plasma, |33j finds that to leading order the 
scalar hydrodynamic QNM is described by 



ujL\, , = ^ V -i— + 0{^], (5.12) 

\hydro ^ g \^^2^y ' V ; 

and this fixes the black point in Fig. |4]when r+/L — ?■ co. This analytical result was already 
compared against numerical data for large radius in [33]. Much like in the vector case, in 
the main sequence {p > 1), microscopic modes with fixed r^/L scale linearly with r^/L. To 
illustrate this property we connect with auxiliary dashed black straight line the two set of modes 
with r_(_/L = 0.2, 0.5. The green dots have exactly the values taken from Table 1 and 2 of [33] 
(we have added the r+ = 0.5 green points in the p > 1 curves). 

In Fig. [5] we give details for the hydrodynamic scalar QNM of the i = 2 harmonic. We 
plot the real (imaginary) part of the frequency as a function of r^/L in the window of values 
where the frequency varies the most, namely r^/L < 7. The red point is the scalar normal 
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Figure 5: Hydro dynamic scalar QNM of the £ = 2 harmonic of GAdSBH. The green dots are 
exactly the values taken from taken from Table 1 and 2 of [33j. The blue (magenta) 
dots are obtained solving numerically the Teukolsky (KI) equations. (See text for 
detailed discussion of these plots). 



mode frequency (5.8) for p = 0. This plot complements the data in Table 1 and 2 of |33j 
which is also represented as green dots in these plots (these are the QNM for r^/L = 0.2, 1, 5). 
For larger values of r^/L the real (imaginary) part keeps decreasing (increasing) monotonically 
(for reference, for = lOOL one has loL ~ 1.732 — 0.007 i [33j). In the limit r_|_/L — )• oo it 
approaches the real value, coL = 1.732, dictated by (5.12) which is plotted as a black dashed 



line for reference. In this Fig. [5j the blue points are obtained solving numerically the Teukolsky 
system with BC (4.18), while the magenta points represent the numerical solution of the KI 



system with BC (4.11). It follows from the map constructed in Section |4] that these two curves 
have to overlap. This is indeed the case (the magenta points are drawn larger to be visible in 
the plot). 

The evolution of the scalar QNM spectrum as the scalar hydrodynamic harmonic changes 
is displayed in Fig. [g] where we plot the curve(s) of Fig. [sj and in addition the hydrodynamic 
QNM of the next seven vector harmonics £. More concretely, from bottom to top (on the Left 
panel and on the right side of the Right panel) we have the harmonics i = 2, 3, • • • ,10. On 
the Left panel we normalize the real frequency to the respective p = normal mode frequency 



of AdS (5.8). We see that the real part of the frequency spectrum is in a first approximation 



isospectral at each fixed r^/L. From the Right panel, we conclude that the imaginary part of 
the frequency always has a minimum, unlike the vector case. Like in the vector QNM case, the 
inset plot shows that there is a window of r^/L around the global AdS case {r+/L = 0; red 
dot (0,0)) where the imaginary part of the spectrum is approximately flat increases, and this 
window increases as the harmonic i grows. The green dots are exactly the values taken from 
Table 1,2 and 3 of |33]. 



As a final remark, note that as discussed below (4.14), the Teukolsky formalism describes 
only the harmonics with i > \s\ = 2. So, it misses the Kodama-Ishibashi vector mode with 
£v = 1 and scalar mode with ig = 0- The QNM spectrum of these KI modes is very special 
because it only contains a zero-mode, i.e. a mode with zero frequency. The scalar zero-mode 
produces a shift in the mass of the solution, while the vector zero-mode generates angular 
momentum (thus connecting perturbatively global AdS-Schwarzschild to Kerr- AdS) |14| I54j . 
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Figure 6: Hydrodynamic scalar QNMs of the first 8 harmonics of the GAdSBH. Viewing from 
the right side of the plots, from bottom to top we have: £ = 2,3,, 10. The red dots 



give the scalar normal mode frequencies (5.8) with p = 0. The green dots are exactly 
taken from Table 1-3 of [SH]. 



The Teukolsky formulation is blind to these modes that generate deformations in the conserved 
charges [271133]. 
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A An overview of perturbations in Kerr(-AdS) 

The Teukolsky solutions and 5^4, the Starobinski- Teukolsky identities, and the Hertz map 
that constructs the associated metric perturbations provide the complete information about the 
most general metric perturbation of the Kerr(-AdS) black hole [H [9l [10] (the only exception 
being the "£ = 0, 1" modes that simply add mass or angular momentum to the background |55j : 
onwards we omit these exceptional modes from our discussion). 

In this appendix we provide a brief historical overview of the studies that culminated with 
the above conclusion. In this discussion we assume the background to be Petrov type D and 
we highlight some facts that are sometimes not duly appreciated. 

In the Newman-Penrose formalism, the fundamental gravitational variables are the 4 com- 



ponents of the NP tetrad basis = {£,n,m,m} , the 12 complex NP spin coeficients (2.2) 
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and the 5 complex NP Weyl scalars {^Oi " " " ^4}- These variables are governed by a set of three 
systems of equations, namely, the Bianchi identities, the Ricci identities, and the commutation 
relations for the basis vectors. The metric is determined once we fix the tetrad basis by 

gf,u = -2£(^n^) + 2m(^mj,) . (A.l) 

The most general perturbation of this system requires determining 10+24+16=50 real func- 
tions to specify the perturbations of the 5 complex Weyl scalars, 12 complex spin coeficients 
and the 16 matrix components that describe the deformations of the tetrad via ^e^ = A^e^,. 

In the Kerr(-AdS) background, this general perturbation system divides into two sectors 

mm- 

I) J^'o, S^^i, 5^3, (5^'4, 6k, 6a, 6\, 6u; (A.2) 
//) 6"^2, 6a, 6(3, 6e, 6^, 6tt, 6p, 6fi, 6t, 6£, 6n, 6m, 6rn. (A. 3) 

The first family describes perturbations of those variables that vanish in the Kerr(-AdS) back- 
ground because it is a Petrov type D geometry. The second involves all other quantities that 
are not required to vanish in such a Petrov background]^ A remarkable property is that these 
two sectors of perturbations "almost decouple" in the sense that we can solve the perturba- 
tion sector /) without solving the NP equations involving the perturbations of sector //). The 
solutions of sector /) are however a prerequisite to then search for the solutions of sector II) 

13 do]. 

Teukolsky [8J, and later Chandrasekhar [S] following an independent computation, found 
the solutions of the perturbation sector I). One just needs to solve the Teukolsky master 
equation that gives the solution for 6'^q, say. The Starobinski- Teukolsky identities fix the 
relative normalization between these two variables [36^ [371 [38l [9l [10]. These scalars 6'^q and 
(5^'4 are gauge invariant, i.e. invariant both under infinitesimal coordinate transformations and 
infinitesimal changes of NP basis. We can then set 6^1 = = 6^^ by an infinitesimal rotation 
of the tetrad basis. Finally, the perturbations of the spin coefficients 6k, 6a, 6\, 6v are obtained 
by applying differential operators to 6'^q and (5^4. So the information on the linear perturbation 
of our system /) is encoded in the gauge invariant variables 6^^^ and ^^'4. (This is not the full 
story concerning perturbation sector I); we will come back and complete it in the end of this 
Appendix) . 

The most general perturbation problem is however not yet solved since we still need to 



find the solutions for the perturbation sector II) in (A.3). In a tour de force computation that 



requires starting with sector /) solutions, Chandrasekhar did a direct and complete integration of 



the remaining linearized NP equations to find the sector II) solutions (A.3) [lOJ. Remarkably, in 
the end of the day, all sector II) perturbations are determined also only as a function oi6^o, 6^4. 
This justifies the statement that the Teukolsky master equations and the Starobinski- Teukolsky 
identities encode the complete information about general perturbations of the Kerr(-AdS) black 
hole. Note in particular that with the knowledge of the basis vector perturbations we can also 



easily construct the perturbations of the metric components through the variation of (A.l ) [lOj. 

An astonishing twist in this story is that we do not need the major effort of integrating the 
full system of NP equations to get what is often the most desired result, namely the metric 
perturbations hab- This was first realized by Cohen and Kegeles |24t I26| and Chrzanowski [25\ 
who have assumed some ad hoc, but smart guessed, hypothesis to build the Hertz map. At this 
point in time, this map was a prescription to reconstruct the most general perturbations of the 
metric tensor only from the knowledge of the Teukolsky master solutions 6^q, 6^4, i.e. without 



^This includes the spin coefficient e that would vanish if we worked with an appropriate tetrad basis, but not 
otherwise. In particular it is not required to vanish by the algebraically special character of the spacetime. 
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requiring information on the variables (A. 3) (in a Petrov type D background). An elegantly 
simple proof of the Hertz map construction was finally provided by Wald in [27J. It promotes 
the Hertz construction from a prescription into a formal map. Keypoints in Wald's proof are: 
1) the existence of the decoupled Teukolsky master equations ( 2.10 )-( 2.11 ) for 5^o, ^^4; 2) 
the fact that the Teukolsky operator for is the adjoint of the one for 5^4, which on the 
other hand is in the end of the day responsible for the existence of the Starobinsky- Teukolsky 



differential identities (2.19)-(2.20); 3) the fact that the equations (2.24) defining the Hertz map 
are the adjoint of the original Teukolsky master equations (2.10)-(2.11 ). Ultimately all these 



properties are due to the algebraically special character of the background. 

The upshot of Wald's proof of the Hertz map prescription of [23]- [26] is that we can use this 



map (2.27)-(2.28) to obtain the complete, most general, metric perturbation of the Kerr(-AdS) 
black hole (with the exception of the modes that change the mass and angular momentum), 
without needing to integrate the extra NP equations that would be necessary to find the solutions 



(A.3). It is worth to look back and appreciate this result: d priori we had to find a total of 50 
variables. However, in the end of the day we just need to solve the Teukolsky master equation 
for d^Q (the equation for <5^'4 is its adjoint and their relative normalization set by the ST 
identities) to get through the Hertz map the most general metric perturbation. An important 
corollary of this result is that the global AdS boundary conditions we find in Section [3] apply 
to generic perturbations of the Kerr- AdS black hole (with i > 2). 

We now return to an issue that was left without its complete discussion. As said above, 
once we have completed the analysis of O |9] we find the solution of perturbation sector I). 
Strickly speaking there is however some residual incompleteness in this solution since [9l [10] : 
1) at this point we just know the absolute value but not the real and imaginary parts of the 
angular St arobinski- Teukolsky constant Cst, and 2) there is still an unknown numerical factor 
needed to fully determine the perturbations 5X, and dv. Both these gaps in our knowledge are 
filled once we solve perturbation sector II) via an integrability condition [9l [10] (also reviewed 
in sections 82 to 95 of chapter 9 of the textbook [H]). To our knowledge, the analogous 
computation of [TO] that determines this information in the Kerr- AdS case was never done, and 
it would be interesting to undergo this task. However, to determine the asymptotically global 
AdS BCs of Section [3] we do not need this knowledge at all. Moreover, we do not need the 
explicit expression for Cst to construct the map between the Kodama-Ishibashi and the a = 
Teukolsky formalisms of Section [4| Nevertheless, and for completeness, in the main text we 
conjectured the expression for Cst to be the solution of |CstP as given in (2.23) that reduces to 
the asymptotically flat expression of [10] when L — )• oo. This expression is written in (2.23) or 



(4.17) when a — )• 0. This is a reasonable expectation but it would nevertheless be important to 
confirm this expression with a computation similar to the one done in [10]. This would close 



the Kerr-AdS linear gravitational perturbation programme. 

It is believed that Einstein's equation is not a special system of coupled PDEs. On the other 
hand, when these equations are linearized around a Petrov type D background and written in 
the Newman-Penrose formalism it is astonishing to find how special the linearized PDE system 



IS. 
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